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Abstract 

We derive a condition that is necessary and sufficient for the instability of inviscid, two- 
dimensional, plane parallel, shear flow with equilibrium velocity profiles that are monotonic, 
real analytic, functions of the cross stream coordinate. The analysis, which is based upon the 
Nyquist method, includes a means for delineating the possible kinds of bifurcations that involve 
the presence of the continuous spectrum, including those that occur at nonzero wavenumber. 
Several examples are given. 
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1 Introduction 

The linear stability of inviscid, incompressible, two-dimensional, plane parallel, shear flow was 
considered over a century ago by Rayleigh, Kelvin, and others. A principal result on the subject is 
Rayleigh's celebrated inflection point theorem [1], which states that for an equilibrium flow to be 
unstable, the equilibrium velocity profile must contain an inflection point. That is, if the velocity 
profile is given by U (y), where y is the cross-stream coordinate, then there must be a point, y = yj, 
for which U"(yi) = 0. Much later, in 1950, Fjortoft [2] generalized the theorem by showing that, 
moreover, if there is one inflection point, then U"'(yi)/U'(yi) < is required for instability (see 
[3] for further extensions). Both Rayleigh's Theorem and Fjortoft's subsequent generalization are 
necessary conditions for instability, but they are not sufficient. That is, even though an equilibrium 
profile may contain a vorticity minimum, it is not necessarily unstable. The point of this paper is to 
derive, for a large class of equilibrium velocity profiles, a condition that is necessary and sufficient 
for instability. 

1.1 Overview 

The procedure we use to derive the instability condition is inspired by techniques developed for 
the Vlasov equation. For that problem one can find a condition that is necessary and sufficient 
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for instability by using the Nyquist method, a method that leads to what is known as the Penrose 
criterion [4] in plasma physics. The Penrose criterion follows fairly straightforwardly for the linear 
Vlasov problem because the discrete eigenvalues satisfy an explicit dispersion relation. However, in 
the context of the Euler equation governing the shearing fluid, the relevant eigenvalue problem leads 
to Rayleigh's equation. This equation cannot be manipulated into an explicit dispersion relation, 
and a sufficient condition for instability of velocity profiles of a general form has not previously 
been given. We note, though, that [5] and [6] present a sufficient condition for long wave instability, 
and in [5] a necessary and sufficient condition for instability of profiles with a single inflection point 
(where bifurcations are restricted to occur through zero wavenumber) was obtained. 

Here we derive a condition that is necessary and sufficient for instability for a class of velocity 
profiles, U(y), where y is the cross stream coordinate. Specifically, we consider profiles that are 
monotonic functions when y G [—1,1], which we refer to as the "flow domain," and which are real 
analytic. This latter limitation means that U(y) has a convergent Taylor series on [—1, 1] and thus 
possesses an analytic continuation into the complex y-plane. Hence, there exists a neighborhood 
(an open set of the complex plane) that contains the interval [—1,1] in which U(y) is analytic and in 
which U'(y) / 0. We assume this neighborhood is as large as is needed in subsequent calculations. 
Both the monotonicity and analyticity limitations can be generalized, but we will not attempt this 
here. 



1.2 Summary of result 

In the remainder of this section we summarize our main result, which amounts to a prescription for 
obtaining the condition that is necessary and sufficient for instability: given the solution, ip(y,c r ), 
to the Fredholm integral equation, 

V>G/,c r ) = g(y,y c ) + £ ^^^U'^y'^iy',^', (1) 

where c r = U{y c ) and Q(y, y') is a Green function (given in Eq. 8 below) containing the streamwise 
wavenumber k and embodying the boundary conditions, we construct the "Nyquist function," 

J-i U{y)-Cr U\y c ) 

where V denotes the Cauchy principal value. We then plot e(c r ) on the (e r , e^)— plane for c r 
running from U{— 1) to U(l), or, equivalently, for y c along the interval [—1,1]. The profile U(y) 
is exponentially unstable for that value of k used in Eq. (1) if and only if the resulting path loops 
around the origin, e r = €j = 0. In fact, the path can only cross the e r — axis at the inflection points, 
yi, of the velocity profile where U"(yj) = 0. This leads to the following necessary and sufficient 
condition for instability: e r (cj) < with cj = U(yj), or equivalently 



'«^>l, (3) 
-l U(y)-ci 



for one of the y/'s. 

In Section 7, among other examples, we treat the profile U(y) = tanh(/%) and arrive at the 
"Nyquist plots" shown in Fig. 1. For this equilibrium profile with (3 = 2 the plots of Fig. 1(a) 
loop around the origin and e r (c[) < (where cj = 0) provided k < 1.832. Hence for (3 = 2 the 
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Figure 1: Nyquist plots for the single inflection point profile, U(y) = tanh/3y. (a) Four plots for 
(3 = 2 and k = 0, 1, and 2, and the critical value for the onset of instability, k = k c ~ 1.832 (dashed 
curve), (b) Four plots for k = and (3 = 1, 1.5, and 2, and the critical value for the onset of 
instability, j3 = (3 C ~ 1.2 (dashed curve). 

profile is unstable over the band of wavenumbers < k < 1.832. Figure 1(b) illustrates the onset 
of instability, which occurs through k = 0, as the parameter (3 is changed. 

In order to prove the result outlined above we need various mathematical results concerning 
Rayleigh's equation. These are derived in Sections 2-5. Application to specific equilibrium velocity 
profiles is presented in Sections 6 and 7. Finally, in Section 8 we summarize, place our work in 
context, and discuss future developments. 

2 Review 

2.1 Formulation 

For a two-dimensional, inviscid fluid contained within the channel, x <G (— oo, oo) and y G [—1, 1], 
an equilibrium state is given by any flow, (U, 0), where U = U(y) is the equilibrium velocity profile. 
Infinitesimal disturbances to such an equilibrium are described by the equation [1], 

d t Co + U(y)d x u> = U"(y)dJ, (4) 

where Q(x,y,t) is the vorticity perturbation, ijj(x,y,t) is the associated streamfunction, which are 
related by 

u = V 2 Vs (5) 

and the boundary conditions are ip(x, ±1, t) = 0. Throughout this paper we will consider solutions 
of the form Co = u(k, y, t) exp(ikx) and ip = i(j(k, y, t) exp(ikx), in which case (4) becomes 

d t uj + ikU{y)u = ikU"{y)^ , (6) 
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and the inverse of (5) takes the form, 



Mv,t) = J^g(y,y'My',t)dy', (7) 

r( ,s J — sinhA;(l — y) sinh k(l + y')/fcsinh2/c if y > y' , . 

^ V,V ' = \-sinhfc(l + y)siiihfc(l-^)/lfesinh21fc if y < y' . ^ 

Here and henceforth we suppress the /c-dependence in the arguments of all functions. 
2.2 Laplace transform 

We need some results associated with the Laplace transform approach to the solution of (6) as an 
initial value problem; so we give a brief review of this approach here. More detailed accounts arc 
given in [7], [8], and with considerably more rigor in [9]. 
Recall, the Laplace transform pair is 

roo 1 r 

*(y,p) = J o e- pt i>(y,t)dt, ^(y,t) = — J J* *(y,p) dp , (9) 

where C is the Bromwich contour that runs parallel to the imaginary p-axis and to the right of all 
singularities of the integrand. Using the above, (6) becomes 

( d 2 2 ikU" \ u(y,0) 

where u(y, 0), the initial value of the vorticity, satisfies u(y, 0) = (d 2 /dy 2 — k 2 )ip(y, 0). The solution 
to (8) can be written formally as 

•<*.')= £*M»)j$g$fi«. (») 

where the Green function, Ti(y,y';p), is given by 

*{y,y,p) \-M, 1 (y>, p )y 2 (y, p )/w(* 1 ,y 2 ) if y'< y <i, {LZ) 

with the Wronskian 

W(*i,tf 2 ) :=^i(y,p)* / 2 (y,p)-^' 1 (y,p)* 2 (y,p). (13) 

Here and ^ 2 are any two solutions to (10) with the right-hand side set to zero and that satisfy the 
boundary conditions, l,p) = ^ 2 (l,p) = 0. Since (10) does not contain a term involving the 
first derivative with respect to y, it is an elementary result (see [10], Chap. V) that the Wronskian 
is independent of y. Arguments pertaining to the inverse Laplace transform of (11) imply that the 
dispersion relation for discrete normal modes is given by 

W(*!,* 2 ) = 0, (14) 

with Re(p) > 0. 



where 
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2.3 Normal modes 

An alternative approach to Laplace transformation follows when we search for normal modes at the 
outset by assuming that u(y,t) = u(y) exp(-ikct) and ip(y,t) = ip(y) exp(— ifcct). Inserting these 
expressions into (6) yields Rayleigh's equation, 

(U-c)u = (U - c)(4>" - k 2 i>) = U"i>. (15) 

which is an equation for the eigenfunction ip{y) corresponding to the eigenvalue c. In the next 
section we describe some properties of this equation, for both real and imaginary values of c, that 
we will need later. 

3 Rayleigh's equation 

3.1 Singular points and solutions in the complex plane 

We will need several facts about the solutions to Rayleigh's equation (15). Firstly, the point y c , for 
which U(y c ) = c, is a singular point of (15). If c is real then this point is commonly referred to as a 
"critical layer." Since U (y) is a monotonic function for y £ [—1,1] there is at most one such critical 
layer in the flow domain, and this occurs when c is in the range of U(y). Because of monotonicity 
and real analyticity of U for y in the flow domain, there exists a neighborhood of [—1,1] in the 
complex plane in which the only singular point of (15) is the critical layer. We add that, about any 
point of the flow domain, there exists a neighborhood in which the complex variables y c and c are 
in one-to-one correspondence. Without loss of generality, we will take U'(y) > when y € [—1,1] 
for definiteness. 

Rayleigh's equation has a singular point at the critical layer, but this is not the only possible 
singular point. First, there may be other points in the complex y plane for which U(y) = c. 
Second, although no others exist when U is continued into a neighborhood of [—1, 1], it is likely 
that additional singular points occur in U when that function is further analytically continued 
into the complex plane. These singularities can then show up in Rayleigh's equation (though it is 
worth noting that because U occurs in Rayleigh's equation only in the coefficient U"/(U — c), for 
meromorphic U all singular points in the finite complex plane are regular; this can easily be shown 
by Laurent expansion). However, here we will only need to consider the critical layer. 

From the elementary theory of ordinary differential equations in the complex plane (see e.g. 
[11], Chap. X) it is known that there are two analytic solutions for y in a neighborhood of any 
ordinary point. Moreover, in general, branch points occur in the solution at the locations of regular 
singular points. If we fix the two parameters, k and c, then the locations of the singular points are 
fixed for a given equilibrium profile. 

In a neighborhood of the critical point in the complex y-plane it is straightforward to obtain 
two Frobenius series solutions of the following form: 

ip g (y,c) = [U(y)-c](p 1 (y,c) (16) 

and 

i>b{y,c) = [U(y) -c]log|E%) -c\ip 2 (y,c) +ip 3 {y,c), (17) 
where (f\, ip 2 , and ips are analytic for y in this neighborhood and for c in the finite complex c-plane. 
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The Frobenius solutions can be used to construct neutral eigenmodes (with c real). We will go 
about this construction shortly; however, Rayleigh's equation may also have complex eigenmode 
solutions, for which c is complex, and we consider these first. 

3.2 Complex eigenmodes 

When c is complex, there is no singular point in Rayleigh's equation for y G [—1,1], and we 
may then eliminate to between (15) and (7) to obtain an integral equation for the streamfunction 
(eigenfunction) of the complex eigenmodes: 



This has the form of a Fredholm equation, and if there is a solution for c = cd, that solution is 
known to be unique (e.g. [12]). Moreover, the conjugate, c = c^,, is also a solution. These complex 
conjugates comprise a pair of discrete eigenvalues, of which there can only be a finite number. In 
fact, at the end of this paper, we essentially give an upper bound on their number. 

The complex pairs that satisfy (18) are equivalent to the homogeneous solutions of (10). Alter- 
natively, the values c = cd and c = c* D are the zeros of the Wronskian (14) of the Laplace transform 
theory. 

According to the Rayleigh-Fj0rtoft Theorem, instabilities occur only if the velocity profile con- 
tains a minimum in vorticity (since we assume U' > 0, the vorticity is everywhere negative). Thus, 
because Im c > signifies instability, the theorem implies that the complex, discrete eigenmodes 
can only exist if the vorticity has a minimum. 

3.3 Neutral discrete eigenmodes 

The solution of Rayleigh's problem for the neutral eigenmodes is not so clear cut. In fact, if c is 
real, then there can be no nontrivial, regular solutions for the streamfunction, with only a single 
exceptional type of mode. This follows from two facts. First, if c < U(—l) = U\ or c > U(l) = U2, 
then Rayleigh's equation with the assumed boundary conditions has only a trivial solution. Thus 
neutral solutions must have critical layers. Second, if c lies within the flow domain, [U\, U2], then 
it can be shown that ip g cannot satisfy the boundary conditions (e.g. [13]). Therefore, all neutral 
eigenmodes must contain the "bad" solution ipb(y, c), in which case ip(y, c) must have an undefined 
derivative. This failure to construct real analytic neutral eigenfunctions is connected to the presence 
of a continuous spectrum that occurs for wavespeeds lying inside the flow domain: c G [U(—l), U(l)] 



The exceptional class of modes are those for which the critical layer lines up with an inflexion 
point of the equilibrium flow profile. In this case, U"(y c ) = 0, and the singular point is removed 
from Rayleigh's equation. The Frobenius solution, ^(y,c), is then an analytic function. The 
exceptional modes also satisfy equation (18), but now c = U(y c ) = U(yi), where yi is the location 
of an inflection point. We call these eigensolutions, "inflection-point modes." These are discrete 
eigenmodes embedded in the continuous spectrum. 

One important feature of the inflection-point modes is that they must be the limits of the 
complex pairs as Im(c) — > 0. In fact, they are the basis of the Tollmien-Lin perturbation theory, 
which builds the nearby complex solutions from the inflection-point mode. This perturbation 
expansion is normally performed with an underlying assumption that the limit exists. Though this 




(18) 



[7]- 
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is largely a technicality, the assumption is strictly only valid when U(y) is an analytic function on 
[—1,1] (this is one of the results that comes directly out of the complex analysis described here). 



3.4 Singular eigensolutions 

Neutral eigenmodes cannot, therefore, be real analytic at the critical layer if U"(y c ) ^ 0. Conse- 
quently, because the derivative of the streamfunction is not defined for these modes, we can only 
impose continuity on such solutions to Rayleigh's equation. Thus tp(y, c) may have an arbitrary 
jump in its first derivative at y = y c . Any family of such solutions can be taken to be a set of 
singular eigenfunctions of the continuous spectrum. However, there is one family that is especially 
useful, and we now construct these special singular eigenmodes. 

The singular eigenfunctions are generalized function solutions to Rayleigh's equation (see [14] 
and [15]): 

"(l/, c) = V^- c + A(c) S[U(y) - c], (19) 

where V indicates the Cauchy Principal Value, S(U — c) is Dirac's delta function, and A is yet to 
be determined. With the interpretation of the singularity in (19) by means of the Cauchy principal 
value, we may define the jump in the derivative of ip by X/U'(y c ). This quantity is currently 
arbitrary and in particular can be taken to be any function of c. Our special singular eigenmodes 
arise from a judicious choice for A. 

If we integrate (19) across the channel, we obtain 



or 



where U' c = U'(y c ) and S(c) is a parameter that is, in general, a function of the eigenvalue. In fact, 
since u is a generalized function, the amplitude of the eigenmode can only be fixed on multiplying 
by a suitable test function and integrating. In Eq. (21), the test function is simply unity, and so S 
can be regarded as the eigenfunction amplitude. Moreover, since this is a linear problem, we may 
choose the dependence of E as we wish; this then amounts to the selection of A. 
On using (7), (19), and (22), we find: 

lKv,c) = ±9<y,Vc)+V l\eM) UH uv)^? ^ (23) 

= H(c)0(y, y c ) + £ g(y, ^^' yc V ( 1 /)V>(y / , c)dy' . (24) 

Although (23) is at first sight a singular integral equation, (24) is a regular Fredholm equation of 
the second kind {e.g. [12]) that is straightforward to solve. 

Provided the Fredholm equation has no homogeneous solution, (24) has a unique particular 
solution. If (24) does have a homogeneous solution, particular solutions are, in general, unbounded. 
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The important point is that a family of well-behaved singular eigenmodes is needed to represent 
the continuous spectrum. If there are no homogeneous solutions, any choice for S(c) will suffice for 
this purpose. But if there are homogeneous solutions, a more specialized choice must be made. One 
possible selection is E(c) = V(c), where T>(c) is the Fredholm determinant. With this selection, the 
inhomogeneous term automatically vanishes if there is a homogeneous solution. This ensures that 
the solution to the Fredholm problem is always bounded. Hence, even if there are homogeneous 
solutions, we can always find a set of sensible singular eigenmodes. (Note that V(y c ) is determined 
by the kernel in (24), that is Q(y,y') and U(y), and is independent of the solution, ip(y,c), and so 
there is freedom for this selection for H(c).) 

For Rayleigh's problem with the profiles we have considered, we have not found any homoge- 
neous solutions to (24). So S(c) = 1 is a convenient choice in any practical application. Then, 
from (19) and (22), we compute the singular eigenfunction and the jump in the streamfunction's 
derivative. 

Finally, note that the amplitude of the solution at the critical layer, ip(y c ,c), cannot vanish, 
since the Frobenius solution (17) satisfies tpf ) (y c ,c) ^ 0. This is an important property of the 
singular eigenfunctions that will be used later. 

4 The dispersion relation 

In this section we construct an expression for the dispersion relation. The form of the dispersion 
relation is designed to facilitate the subsequent Nyquist analysis, and is one that is useful for 
relating discrete eigenmodes to continuum eigenmodes. 

4.1 Discontinuity and the dispersion relation 

We begin by considering an arbitrary point y* that is not coincident with the critical layer; i.e. y* 
is any ordinary point of the differential equation. For convenience we choose 6 [—1,1] and we 
assume that c is fixed in the upper half c-plane. 

Now we construct two solutions, \t<(y, c) and \I/>(y, c), that are defined for y £ [—1, 1] on the 
complex c-plane. The first, ^f<(y, c), is defined by series expansion about the left boundary point 
y = — 1 and satisfies the boundary condition \I/<(— l,c) = 0. This is a one-parameter family of 
solutions, where the parameter can be taken to be a multiplicative constant; i.e. we can write 
^<(?/; c ) = C < $ < (y,c), where $<(y, c) is a parameter-independent solution of Rayleigh's equation 
that satisfies $<(— l,c) = and $><(— l,c) = 1. Similarly, a second one-parameter family of 
solutions is constructed by series expansion about the point y = 1. We denote this second solution 
by *>(y,c) = C><l>>(y, c), where by construction $>(l,c) = and $>(l,c) = 1. We emphasize 
that at fixed y £ [—1, 1], both ( ! > <(y, c) and $>(y, c) are analytic in c for all c ^ U(y). 

The next step in our construction is to choose the constants C < and C> so that ^f < and ^ > are 
continuous at the point y*; i.e. we set *<(y*,c) = \I/>(y*,c). This requirement leaves a remaining 
constant that is an overall scaling factor; the two solutions can be represented as follows: 

*<(y,c,y*) = C(c,y*)$>(y*,c)$<(y,c) 

c,y*) = C(c,y*)$<(y*,c)$>(y,c), (25) 

where the sole remaining constant is C(c, y*), in which we have included dependence upon c and 



8 



y*. Clearly, we are free to choose C(c, y*) arbitrarily. We give a prescription for denning C(c, y*) 
shortly 

No matter which (nonzero) value is chosen for C(c, y*), the derivatives of the functions '!'< 
and ^> will not, for general values of c, match at the point y*. However, in the event that they 
do match, the value of c is an eigenvalue and the functions and ^> define an eigenfunction. 
In fact, when this is the case, ^> < and x I / > are analytic continuations of each other. This follows 
because y* is an ordinary point and both ^< and ^> are solutions of the Cauchy problem for 
Rayleigh's equation with identical specification of their values and derivatives at y*. Moreover, 
when the derivatives of the solutions match, the point y* is immaterial and the dependence upon 
y* drops out of the incipient eigenfunction. 

In light of the above, the jump in the derivatives of and \I>> at y* is equivalent to a dispersion 
relation: its vanishing determines c = c(k). (Recall that the k dependence has been suppressed in 
the expressions above.) Note that if C(c, y*) is chosen to be an analytic function of c in the upper 
half plane, then ^<(y^,c) and *I' > (y i(i ,c) are also analytic functions of c in the upper half plane, 
which follows from the assumption U(y*) ^ c. Hence, the dispersion relation is an analytic function 
for c in the upper half plane. Also, note that the dispersion relation does not depend upon the 
point y*, since matching of the derivatives at any point will give an eigenfunction. 



4.2 Relationship between e and W 

We now obtain a convenient expression for the dispersion relation. 

In general, for y within the flow domain, we write the streamfunction in the form, 

*(y,c,y*) := H(y - y^^y^cy*) + H(y* - y)* <(y,c,y*) (26) 

and H is the Heaviside function. The vorticity, by which we mean w(y, c, y*) = — k 2 ^, then has 
a representation, 

w(y, c, y*) = Q(y, c, y*) + e(c, y*) <5(y - y*) . (27) 
where e(c, y*) measures the jump in the derivatives, 

e(c,y*) := *>(y*, c, y*) - #<(y*,c,y*) , (28) 

and 

rv w r „ ^ - J U"(y)*<(y,c,y*)/[U(y) - c] if y > y* 

S ^' C ' y * j - \C/"(y)* > (y,c,y,)/[C/(y)-c] if y < y* . ^) 
By inserting (25) into (28) we see that 

e(c,y*) = C(c,y*) [$>(y*,c)$' < (y*,c) - $<(y*, c)$>(y*, c)] 

= :C(c,y*)W(c), (30) 

where VF, the Wronskian for Rayleigh's equation, is independent of y*. Note that this quantity can 
be identified with the Wronskian of (13) by substituting p = —ike into the latter. 

We next normalize the solution by integrating (27) over the flow domain and setting the result 
equal to the c— dependent parameter, S(c, y*): 

J i ^{y,c,y*)dy = E(c,y*) , (31) 
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which yields 

E(, ri = ifet ) + r gWtSt) %+ fy^' * (32) 
y_i U{y)-c J yt U{y)-c 

By using Rayleigh's equation, we then find 

E = -*'<(-!, c, + *' > (l,c,y*) - fc 2 T* *<(y,c,y,)dj/ - fc 2 *>(y, c, y*) dy . (33) 
Finally, by inserting (25) into (33) we may express C(c, y*) in terms of S(c, y*) (or vice versa): 
C{c,y*) = -E(c,y*)\k 2 [ V $>(y*, c)$<(y, c) dy + /c 2 / $<(y*, c)$>(y, c) dy 



-i 



+ $>(*/*> c) -$<(y*,c)J . (34) 

Because <3?< and <£> are analytic in c in the upper half plane for fixed y, the quantity in the 
denominator of the right-hand side of (34) is also analytic. Thus, if we were to select C(c, y*) to 
be analytic in the upper half c— plane, then S(c, y*) would also be. However, the converse is not 
quite true: if S(c, y*) is selected to be analytic in the upper half plane, then C(c, y*) is also analytic 
except, perhaps, for poles at the zeros of the denominator. We will interpret these zeros shortly. 

4.3 Rayleigh Green function and singular eigenfunctions 

The next piece of the puzzle is to relate the dispersion relation to the singular eigenfunctions. 
We may rewrite (27) in the form, 

*"-fc 2 *-7^ = 6«5(y-y*). (35) 

Thus, ^ is related to the Green function for Rayleigh's equation. This equation is a bit subtle, 
since e (according to (32)) is in fact a property of the solution and so the equation is self-referential. 
However, if y* / y c and e / 0, then \I//e satisfies the normal equation for the Green function. 

We may turn equation (35) into one of integral form on using the Green function of the Lapla- 
cian: 

*( y ,c,y*) = g(y,yMc,y*) + £ G(y,y') U '' {v ^'l C ' c V * ] d y' ■ ( 36 ) 

Moreover, on using the relation (32), this can be written as 

*(y, c, y*) = E(c, y*)G(y, y*) + J *(y', c, y*) /C(y, y', y*, c) dy' , (37) 

where 

fc(y,y ,V*,c) :=U (y ) _ . (38) 

Hence the Green function can be constructed by solving another Fredholm integral equation. 



10 



In fact, the two Fredholm problems are closely related: if we take the limit c — > c r + iO = 
U(y c ) + iO with y* = y c , then we recover the Fredholm equation for the singular eigenfunctions 
from (37)-(38) with ^(y,c r + iO,y c ) = ip(y,c r ) and H(c,y c ) = E(c). 

Similarly, we may recover the singular eigenfunction (19) from the Green function equation (27) 
in the same limit: we let c — > c r + iO = £7(y c ) + zO in (32), giving 

~( jl n ^ ^7 j_ n ^ ?7 // (j/)^(y,c r + i0,y c ) 

e(c r + zO,y c ) = =(c r + ?0,y c ) - / — dy. (39) 

However, since \l/(y, c, y*) is analytic in the upper half plane, there exists a generalized form of the 
Plemelj relation (e.g. [16]), which leads to 

, ^ f 1 U"(yU(y,c r ) U"(y c U(y c ,c r ) , , 

=: e r (c r , y c ) + ie;(c r , y c ), (41) 

on using the association of the Fredholm equations to replace Vl/(y, c r + i0,y c ) with ijj(y,c r ). The 
Plemelj relation also implies that 

U"ib(v c ) 

u(y, c r + iO, y c ) = V rJ + e r (c r , y c )5{y - y c ). (42) 

Finally, we identify u>(y,c r + i0,y c ) with to(y,c r ), e(c r ,y c ) with e(c r ), and e r (c r ,y c ) with \{c)/U' c . 
With these associations understood, we will drop the tildes in the next section and the extra 
argument in y c . 

Note that we could take a completely different approach to the problem beginning from the 
Green function. First we select S(c,y*) to be a function that is suitably analytic in c. Then, 
the solution of the integral equation (37) can be shown to have the various analyticity properties 
we have found for ^f(y, c, y*) by using Fredholm theory (the kernel is an analytic function of the 
parameter c in the upper half plane). From there, we build the quantity e(c, y*), which we know 
has zeros if c is an eigenvalue. Moreover, equation (34) determines a function C(c, y*) by which we 
may relate the solution of the integral problem to the solutions <3?< and <&> used above. 

The only fly in the ointment is the possibility of a homogeneous solution to the Fredholm 
problem for some value of c, in which case the particular solution we seek is unbounded. From 
Fredholm theory we know that the solution ^ is not analytic at this value of c, and, in fact, has a 
pole. 

However, this nonanalyticity is evidently an artifact of solving the Fredholm problem, since the 
solutions and are analytic in c. Moreover, since they are related to the solution ^ of (37) 
simply by the factor C(c, y*) as in equation (25), it is clear that the nonanalyticity in \I/(y, c, y*) 
is equivalent to a pole in C(c, y*). But since we may choose E(c, y*) to be analytic, the pole must 
arise from a zero in the denominator of (34). This is the advertized interpretation of the poles in 
C(c, y*); they correspond to the presence of a homogeneous solution to the Fredholm problem (37). 

At this stage, two remarks are in order. We have already noted that there are no homogeneous 
solutions for c in the flow domain in the context of the singular eigenfunction equation. Hence, 
C(c, y*) must be real and contain no poles for these values of c. Secondly, since this nonanalyticity 
is Vl/(y, c, y*) is purely an artifact of solving the Fredholm problem (37), we can take a slightly 
different tack and make a judicous choice for H(c, y*) that avoids the problem. Again, this is just a 
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choice like S = T>(c, y*), where D(c, y*) is the Fredholm determinent of the kernel in (37). Once we 
make this choice, we may derive a solution of the integral equation that is analytic in c in the upper 
half plane, and use it to build the dispersion relation through e(c, y*). However, we will explicitly 
follow the route outlined earlier in this section rather than this Green-function based avenue. 



5 Nyquist method 

In the previous section we discussed the functions e(c) and W(c), and their relationship. If their 
exists a discrete value of the phase velocity, c = cd, such that W(cd) = 0, then we have an 
exponentially growing eigenmode with growth rate, Im^kco)- By construction we know that W(c) 
is analytic for c in the upper half plane, and that it also has a branch cut along the real axis 
between U\ := U(—l) and U2 '■= U(+l). This branch cut arises from that of the natural logarithm 
of the Frobenius solution and the matching procedure (cf. (75) of the Appendix.) Because of these 
properties we begin with W in our analysis below, but in the end we express the final result in 
terms of e, a quantity that is by design reminiscent of the dispersion relation of Vlasov theory. 

The Nyquist method relies upon the argument principle of complex analysis. In the present 
context this principle implies that the integral, 



5^/^* <«> 

2m Jc W(c) v ; 



counts the number of zeros of W in the region enclosed by a closed contour C in the c— plane. We 
choose C to run along the c r axis, with a fixed to an arbitrarily small positive value, and then the 
contour is closed by a large semicircular portion as shown in Fig. 2. As the radius of the semicircle 
goes to infinity, (43) gives the number of exponentially growing eigenvalues. Equivalently, (43) is 
the number of times the path determined by the function W(c) encircles the origin of the (W r , W{) 
plane as c completes a circuit of C. 

It is efficacious to decompose the integral (43) into the following two pieces: 



1 1 w'{c) dc _ J_ f U2+i0 W(c) ± 1 r W'{ 
2ri Jc W(c) C ~2^i J Ul+ io W(c) C+ 2^iJc W\, 



W'(c) 

[ ^dc. (44) 



The first integral is all important, while the second is relatively minor. To understand this recall 
Howard's semicircle theorem [8], which states that the zeros of W must lie within a disk of radius, 
(U\ + U%)/2, centered at the point (U2 — U±)/2 (as illustrated in Fig. 2). Hence, if C encloses the 
semicircle, it contains all of the unstable eigenvales. In fact, wherever W is analytic, we may deform 
the contour C'. Since W is analytic in the upper half plane we may deform C' into any contour 
there that connects W(Ui + iO) to W(U2 + i0) The important point is that since C is deformable 
to any other contour lying outside Howard's semicircle, the path defined by W(c) as c varies along 
C cannot lead to a new enclosure of the origin because this would mean a zero of W outside the 
semicircle. In other words, the count of unstable eigenvales must be independent of the integration 
around C; the only importance of the integral over C is to complete a closed path in the (W r , Wi) 
plane without encircling the origin. In the Appendix we demonstrate this explicitly by extending 
the contour C to infinity and then by analyzing the image of its various pieces in the VF-plane. The 
upshot is that we may ignore the C part of the contour: the change in the argument of W along 
the flow domain is equal to the number of times the function W(c) encircles the origin, which is 
equal to the number of unstable eigenvalues. 
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4 5 



u. 



Figure 2: The contours C, C, and H in the c-plane. The closed contour C runs the entire circuit 
from 1— > 2 — ► . . . — > 8 — > 1, with the portion from 3 — > 1 lifted infinitesimally above the real axis. 
The contour C is C with the piece along the flow domain, 5^7, removed. The contour H (dashed) 
is Howard's semicircle, within which the unstable eigenvalues must lie. 



In Section 4 we showed that e and W are related by e(c) = C(c) W(c). But as c traverses the 
flow domain: 

rU 2 +iO e > rU 2 +iO fU 2 +iO Q' 

\ ~dc= — dc + 



Ui+iQ C 



dc . 



Upon introducing 
(45) becomes 



e e 



i arg 



[e] and W= \W\e iav8lw] , 



(45) 
(46) 



^ In |e| + i arg[e 



U 2 +i0 
Ui+iO 



(\n\W\ +lnC + iarg[VF]) 



U 2 +i0 
U!+iO ' 



(47) 



Now, as remarked at the end of Section 4, C(c) is real and finite along this part of the contour. 
Hence, equating imaginary parts gives Aarg[e] = Aarg[W]. 

Therefore, the number of unstable eigenvalues is determined by the change in the argument of 
the following quantity: 



e(c r + iO) 



1 U"(y)ii(y,c) 



U(y) - U(y, 



dy 



it: 



U"(y c )iP(y c ,c) 
U'(y c ) 



(48) 



(49) 



where c r = U(y c ) (and we have explicitly made the choice S(c) = 1), which is the "Nyquist function" 
of Eq. (2). As c r varies from U\ to U2, or equivalently, y c from —1 to 1, e(c r ) executes a closed path 
in the (e r ,ej) plane. This path begins and ends at (1,0) since ^>(±l,c) = 0. In between, the path 
circulates around and the following considerations help us determine whether or not the origin is 
encircled. 

The path can only cross the e r — axis at the points for which = 0, or ip(y c , c r )U"(y c ) = 0. 
However, as noted at the end of Section 3, ^(yc,c r ) cannot vanish. Hence the crossing points 
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are just the inflection points of the velocity profile. We denote these by yi, with U"(yi) = and 
c = a = U(yi). Thus, the path emerges from the point (1, 0) on the (e r , e^) plane, circulates around 
crossing the e r — axis as many times as there are inflection points, and then terminates at (1,0). 
If the origin is encircled, we have an exponentially growing instability. Hence, by calculating the 
singular eigenfunctions from the Fredholm problem (24) and then constructing e(c r + iO), we can 
determine if there is an unstable eigenvalue. 

In fact, the path can only enclose the origin if there is at least one crossing point to the left of 
the origin. Such a crossing point is guaranteed if 



for at least one of the inflection points. 

Inequality (50) is actually a sufficient condition for instability. We see this as follows: if there 
are an odd number of crossings to the left of the origin, then (50) certainly guarantees an enclosure 
of the origin. But if there are an even number, then one can envision paths that cross to the left 
of the origin, but do not enclose it. In this case, (50) may still hold, but the locus fails to encircle 
the origin. However, the condition still predicts instability because we may yet vary k, and the 
path must change continuously as we change k. From the Fredholm problem, one can observe that 
ip ~ 1/k for large k. This means that ip — > as k — > oo. Thus for large enough k, the path of 
e shrinks about the point (1,0). Since varying k also cannot destroy crossing points, it must be 
the case that there is a range of values of k for which two of the crossing points must straddle the 
origin, and the path encircles it. (An example exhibiting this feature is given in Section 7.3.) 

Thus (50) guarantees an enclosure of the origin for some value of k, and for any number of 
crossings to the left of the origin. In other words, we deduce the following necessary and sufficient 
condition for instability: 

Rayleigh's equation possess an unstable eigenmode if and only if there exists a wavenum- 
ber, k, and at least one inflection point, yj, such that Eq. (50) holds. 

Note that, by "crossing point" we mean strictly only those inflection points for which U'" ^ 0. 
This implies that intersection of the path with the axis leads to a traversal of the axis from one 
half plane to the other. Inflection points for which U'" = can lead to the path touching but not 
crossing the axis, and these intersections will not signify a nearby enclosure of the origin. However, 
an arbitrarily small deformation of U(y) can in this situation lead to the locus encircling the origin. 
Thus profiles with U'" = correspond to marginally stable states (see Section 7). 

Unfortunately, the sufficient condition implied by (50) has a more complicated dependence on k 
than in the related Vlasov problem (the k— dependence in the Vlasov problem is contained purely 
in a multiplicative factor). The dependence here comes from the dependence of ip(y, cj) on k, which 
is not obvious. Nevertheless, if k 2 > k^ := — Min [U"/(U — cj)] for all of the inflexion points yj, 
then we can apply a result of Howard [17] which indicates that there can be no instability. Hence 
for k > k m , the crossing point must lie to the right of the origin, so we can at least bound the range 
of interesting wavenumbers from above. 

Finally, the condition (50) refers only to exponentially growing eigenmodes and, therefore, 
does not encompass all the possibilities for instability. It fails to account for a possible algebraic 
instability of the continuous spectrum. However, in the present context, we can safely ignore such a 
possibility since we know that the continuum is stable for a monotonic velocity profile [9]. Therefore, 




(50) 
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algebraic instability is ruled out. Hence (50) is the encompassing condition for instability. This 
may not be true if the profile is nonmonotonic (see [18]). 



6 Special cases 

In the previous section we described our necessary and sufficient condition for instability; here we 
illustrate our criterion for two special examples: general profiles with k = 0, and an asymptotic 
result. 



6.1 General k = eigenmodes 

For an arbitrary profile with k = 0, Rayleigh's equation simplifies and the eigenmodes can be found 
by directly integrating. It is straightforward to construct the Wronskian, 

W(c) = (U 2 - c)(c - ETi) £ Jjj^- (51) 

We could analyze this dispersion relation independently of the machinery developed in the previous 
sections. However, we consider the k = problem by way of illustration, and so we apply the general 
methodology. 

The Fredholm problem (or the procedure of Section 4) for the streamfunctions of the singular 
eigenfunctions has in this case the analytical solution: 

h( ) = f-(^2-c)(c-i7 1 )(i7 2 -C/ 1 )- 1 ^(C/-c)^ 1 cZy7[C/(y')-c] 2 for y < y c 
™> \ (U 2 -c)(c-U 1 )(U 2 -U 1 )- 1 U&U-c)fidy'/[U(y')-c] 2 for y > y c . { ' 

Consequently, the Nyquist function can be written in the form, 

< i (\\ (^2 ~ c r )(c r - E7j) , f l dy 



or 



(U 2 - c r )(c r -Ui) TT ,{ 1 . 1 . v f 1 U"(y)dy \ 
r[ r) " ' U 2 - C/i c I (U 2 - c r )U> 2 + (c r - U,)U[ + U [U(y) - c r ]U'(yf } [ ^ } 



and 

(U 2 -c r )(c r -U 1 )U" (y c ) 
(U 2 - U^U'iyc) 3 

where U[ = U'(-l) and U' 2 = U'{1). 

Note that, for the k = problem, C(c) = —U' C /(U 2 — U\) for c in the flow domain. Thus C(c) 
is indeed real and finite, as we remarked earlier. 

In the Nyquist construction, the important piece of the contour C leaves the point (1, 0), executes 
some path that is dependent on the details of U(y), and finally returns to (1,0). That is, the path 
is closed, as we remarked earlier. The Nyquist function of (53) is related to that obtained in [5], 
where this special case of k = was considered in a more specific fashion. Note, however, that 
the Nyquist function of [5] differs from (53) by a factor in front of the integral that leads to e(c) 
vanishing at c = U\ and U 2 . If this were not so, our Nyquist plots would not be closed loops and 
one would be forced to consider the neglected piece of the contour, C", in detail. 
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6.2 An asymptotic result 

Our second example is only an approximate result; it concerns velocity profiles of the form, U (y) = 
y + e 2 U(y/e), where e <C 1. This kind of velocity profile represents a linear background profile with 
a superposed, sharply varying, "defect." To leading order in e, the corresponding streamfunction 
is given by 

^(y,c)=g(y,y c ) + 0(e), (56) 
and the Nyquist function has the simple form 

e ( C ) = l-^r^M^ + 0( £ ), (57) 

2k J-oo 7] — c 

or 

, , 1 tanhfc f°° W'Mdr) , 
e r (c r ) = l — — T 7 / — ^^ + 0(e) (58) 



2k J _oo Tj c 

and 



r 



e i (c r ) = - l -^-U"(y c ) + 0(e). (59) 

This asymptotic result is given in [19], following [20]. It is closely related to the corresponding 
Vlasov solution. 

Note that in the two examples, U(y) need not be analytic on [—1, 1]; in fact, the existence of 
two derivatives is sufficient. 



7 Numerical solutions for sample profiles 

We now construct some Nyquist plots numerically for a trio of sample profiles that illustrate different 
features of the inviscid stability problem. 



7.1 Single inflection point 

The profile, 

U(y) = tanh/fy, (60) 

is an example of an equilibrium with a single inflection point. Nyquist plots for various values of k 
and (3 = 2 are drawn in Fig. 1(a). For wavenumbers k < k c with k c ~ 1.832, there is evidentally an 
unstable eigenvalue, since the Nyquist plots enclose the origin over this range of k. Plots for k = 
and different values of (3 are shown in Fig. 1(b). 

In Fig. 1(a), the Nyquist plot passes through the origin without encircling it for the critical 
value k = k c . This wavenumber is the demarcation between stability and instability; that is, it is 
the stability boundary, k c = k c ((3). This stability boundary is displayed in Fig 3(a). As (3 decreases, 
the critical wavenumber decreases until it vanishes at a special value of (3 = (3 m ~ 1.2. For the 
profile with this critical parameter value, the k = Nyquist plot passes through the origin and 
none encircle it. In other words, such a profile is a marginally stable state. 

This feature of the profile is an example of a general result: instability always sets in first at 
k = for profiles with a single inflection point {e.g. [8]). Hence the instability condition is given 
by e r < with e r given by Eq. (54). This is equivalent to the result of [5]. 
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(a) (b) 




Figure 3: Stability boundaries for the three equilibrium profiles of Sec. 7. Figures (a) and (b) depict 
the nucleation of instability through k = 0, while (c) depicts the nucleation of instability through 
k 0. In (b), the stability boundary of the second and third unstable modes is also shown. 
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In some applications, however, one may be interested in flows that are periodic in x (such as 
in annular or spherical geometry, or in numerical simulations). In these cases, there is a minimum 
wavenumber, and k = is neither accessible nor relevant. Hence the k = theory is not applicable 
even for a single inflexion point. In this circumstance one must deal with the general Nyquist 
function and sufficient stability condition described in Section 5. 

7.2 Multiple inflection points 

The profile, 

U(y) = y + —sin-fy, (61) 

is an example in which there are multiple inflection points. In this case, the number of inflection 
points varies with 7. For 7 = 5, 7.5, and 10, there are 3, 5, and 7 inflecton points, respectively. As 
7 increases, the new inflexion points appear in pairs through the boundaries, y = ±1. The Nyquist 
plots at the three parameter values and k = 1 are shown in Fig. 4(a)-(c). For these plots there are 
multiple encirclings of the origin, signifying multiple instabilities. For example, when k = 1 and 
7 = 10, there are three loops around the origin (see Fig. 4(c)). As 7 increases and inflexion points 
appear at the boundaries, the Nyquist curve acquires more loops that appear out of the asymptote 
(1,0). 

At 7 = 10, there are therefore three unstable eigenvalues with k = 1, which is the most this 
profile can support (four of the inflection points lead to crossing of the e r — axis, but these lie to the 
right of the asymptote (1,0)). Note that the antisymmetry of the profile means that the path of 
e(c r + iO) is symmetric under reflection about the e r — axis. 

If we vary k, and calculate e r (cj), where cj with / = 1,2, 3 denotes the three inflection points 
that lead to crossings to the left of the asymptote (1,0), then we obtain the picture shown in 
Fig. 4(d). This shows that the e r (cj) increase monotonically with k. This suggests that if we were 
to vary 7, then all three unstable eigenvalues would appear first at k = 0, which is indeed true as 
can be seen from Fig. 3(b). Thus, the marginal state for this profile is again given by the k = 
theory. 

7.3 Finite wavenumber instabilities 

The third profile, 

U(y) = y + 5y 3 + / tanh 4(y - 1/2) , (62) 

contains one or three inflection points, depending on the value of the parameter /. More specifically, 
as we increase / through about 1.59, we create two inflection points near the point y = 0.6 (see 
Fig. 5). The critical profile for which the two inflection points emerge contains a point with 
U" = U'" = 0. This leads to a Nyquist plot that contains a nontransversal intersection of the path 
with the e r — axis; that is, the plot touches the axis but does not cross it. For larger values of /, 
this degenerate point splits into the two inflection points and the plot crosses the e r — axis twice. 
Nyquist plots beyond this bifurcation are shown in Fig. 6(b). 

This kind of a change to the profile leads to a situation in which there can be two crossings of 
the e r — axis to the left of the origin, but in neighborhood of k = there are no unstable eigenvalues. 
Yet, as we remarked in Section 5, increasing k leads to a deformation of the path such that the 
origin is eventually encircled for some range of k (see Fig. 6(a) and (b)). In Fig. 7, we show the 
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Figure 4: Nyquist plots for the multiple inflection point profile, U(y) = y + (siwyy)/2^, with k = 1. 
(a) The case 7 = 5 has one encircling of the origin and thus one unstable eigenvalue, (b) In the 
case with 7 = 7.5 two more loops have appeared and the curve almost encircles the origin two 
more times. That is, it is close to the threshold of the instability of two further modes, (c) The 
case 7 = 10 has three encirclings of the origin and three unstable eigenvalues, (d) Depiction of the 
monotonic dependence of the crossing values to the left of the asymptote (1,0) as a function of k. 
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Figure 5: Vorticity and vorticity gradient profiles corresponding to the velocity profile U(y) = 
y + 5y 3 + /tanh4(y — 1/2) with equilibrium parameter values / = 1.3, / = 1.65, and / = 2. At 
/ = J c k, 1.59 two inflection points emerge. 
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Figure 6: (a) Nyquist plots for the profile U(y) = y + 5y 3 + /tanh4(y — 1/2) for k = 0, 0.75, and 
1.5, and / = 1.62. (b) Magnification that indicates the onset of instability through k ^ 0. 
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Figure 7: Plots of e r (cj) against k at / = 1.62, where cj, / = 1,2, are the two inflection points 
leading to crossings of the Nyquist plot in Fig. 6 to the left of the asymptote (1,0). 
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variation of e r with k for the two newly created inflection points at / = 1.62; the range of unstable 
wavenumbers is 0.12 < k < 1.31. Again the variation of e r with k is monotonic. 

The onset of the unstable band arises at the value of / for which the two inflection points 
emerge; that is, this profile is the marginally stable state. Moreover, the critical value of k is that 
required to make the nontranversal intersection of the path with the e r — axis occur at the origin. 
Thus the emergence of the two inflection points as we raise / leads to an instability with an onset 
at finite wavenumber, as shown in Fig. 3(c). Thus, in this case, instability is bounded away from 
k = 0, in contrast to the previous examples. (This particular example is somewhat analogous to 
the bump-on-tail instability in the Vlasov problem [4]). 

In both Figs. 4(d) and 7, e r appears to be a monotonically increasing function of k. If this 
feature were generally true, then the deformation of the Nyquist plots as we vary k would certainly 
be simpler to understand. Indeed, this would imply that e r (c/) < evaluated at k = was sufficient 
for instability, independently of the number of inflection points (at least for monotonic, analytic 
profiles). However, we have found no argument implying such a property of e r , and doubt one 
exists. 

8 Closing remarks 

We have presented in this paper a necessary and sufficient condition for the instability of monotonic 
shear flows. Several examples of equilibrium velocity profiles were treated in Sections 6 and 7, 
demonstrating that the condition is of practical utility for finding instability and for understanding 
the kinds of bifurcations that can occur. 

One may ask the question of how the method we have presented, which entails solving a Fred- 
holm integral equation, compares with directly calculating the solutions numerically. Without the 
Nyquist method, one could solve Rayleigh's equation in order to locate unstable modes with a given 
k for c in the complex plane. This amounts to solving a boundary value problem with singular 
points. With this procedure one would repeatedly solve Rayleigh's equation at different values of 
(complex) c and employ some kind of two-dimensional searching technique in order to find the 
eigenvalues. Whilst not especially difficult, this approach is computationally intensive, and in our 
opinion is neither elegant nor insightful. As a computational exercise, the calculation of the Nyquist 
function is substantially less effort than solving the boundary- value problem in this fashion. The 
Nyquist recipe is simpler because numerically it only requires a matrix inversion and only values of 
c on the real axis are involved. Moreover, the sufficient condition of (50) requires one computation 
per inflection point, for each value of k, and no search in the complex plane. 

In addition to these computational advantages, the Nyquist method provides insight into inviscid 
instability. It allows one to determine the number of unstable eigenmodes of a profile, as indicated, 
for example, by the multiple encircling of the origin of the Nyquist plots of Fig. 4. Also, it leads 
us to understand the various kinds of bifurcations to instability of a family of shear flow profiles. 
Put another way, Nyquist theory tells us how we can construct profiles with certain kinds of 
instabilities. For example, with the Nyquist imagery, we know how to build profiles that suffer the 
onset of instability at finite wavenumber, as in Section 7.3. Without this machinery, it would be 
much more difficult even formulating the problem of how to fashion the needed U (y) . 

The bifurcations considered here are ones in which unstable eigenvalues appear out of a con- 
tinuous spectrum. It is important to distinguish this kind of bifurcation from those that occur in 
systems with only discrete spectra, since any stable profile is arbitrarily close (in L p norm) to an 
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unstable one (this is the essence of Gill's result [20] for Couette flow), a feature intimately linked 
with the presence of the continuous spectrum. It would then appear that catagorizing bifurcations 
to instability is problematic. However, the Nyquist plots in the (e r , ej)— plane provide a measure 
of how close a given equilibrium is from being deformed to one at the onset of instability. This 
onset may be through k = (cf. Figs. 1 and 4) or at finite wavenumber (cf. Fig. 6). Indeed, it can 
occur via more degenerate kinds of bifurcations, and these may be visualized straightforwardly with 
Nyquist theory. In fact, it is only the unclear dependence on k that prevents us from immediately 
classifying every instability of a profile with a given number of inflection points (cf. [19]). 

An example of the kind of insight provided by the Nyquist method is given by the following 
general result: for given k, in order to encircle the origin, the path must cross the e r — axis twice. 
That is, in order to create a new instability we need two new inflection points. Thus, if there are 
N inflection points in the profile (and so there are N + 1 crossing of the e r — axis altogether), there 
can be at most (N + l)/2 unstable eigenmodes with that value of k. This result, which builds on 
a theorem of Howard [17], was stated without proof in [21]. 

In finite dimensional Hamiltonian systems bifurcations are regulated by Krein's theorem, which 
states that a necessary condition for the bifurcation to instability is that colliding eigenvalues possess 
energy signatures of opposite sign. It is also know for Hamiltonian systems that the constancy of 
energy can be used to obtain a sufficient but not necessary condition for stability. The results of this 
paper lend insight to and can be interpreted in the Hamiltonian context: the sufficient condition 
for stability based on energy is equivalent to the conditions of Rayleigh and Fj0rtoft, and a version 
of Krein's theorem (and a notion of signature) for bifurcations involving the continuous spectra can 
be related to the Nyquist function [15]. 

We conclude with some remarks about the assumptions we made on the forms of the profile. 
The analytic structure of the Nyquist function heavily relies on the analytic form of the profile in 
the vicinity of [—1,1]. However, irrespective of the form of U(y), we may nevertheless construct 
e(c r ). In fact, the zeros of this function for c r € [—1,1] still point to the existence of regular neutral 
eigenmodes even when U (y) is not analytic. What is no longer clear is how to continue e off the 
real axis. Or equivalently, whether there are discrete complex modes at nearby parameter values. 
Nevertheless there are various hints in this work that suggest that the theory works even if U(y) 
has only two derivatives. 

Finally, nonmonotonicity of U(y) brings in the new complications of multiple critical layers for 
a given wave speed, and unusual properties of the streamfunction for critical layers lying at the 
shearless points where U'(y) = 0. It is for these reasons that the theory does not immediately carry 
over. We are currently working on these and other extensions. 
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Appendix 

In this appendix we show that the change in the argument of the Wronskian is determined entirely 
by the piece of the contour C that is along the flow domain (5 — ► 7 of Fig. 2); the remainder of the 
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contour will be shown to not encircle the origin. 



The big semicircle: 1 — > 2 — > 3. 

On the big semicircle c = i?e* e with 9 € (0,7r). We consider the limit R = |c| — > oo by invoking 
the principle of permanence (see e.g. [22]), which in the present context simply states that the 
limit c — > oo of the solution to Rayleigh's equation at fixed y is equal to the solution of the 
c — > oo limit of Rayleigh's equation. The latter limit is ■0" — /c 2 V> = 0, which has the following 
solution with the appropriate boundary conditions: 3>^?(y) := sinh[fc(y + l)]/k. The principle 
of permanence implies 3><(y,c — > oo) = <E>^?(y), and thus using W(c) = — 3><(l,c) we obtain 
Hindoo W(c) = — sinh(2/c)/fc < 0. Therefore, the piece of the contour 1 — > 2 — > 3 maps into the 
W-plane as a single point on the negative real axis (this is the asymptote (1,0) on the e— plane). 

As a check consider the limit k — ► 0. Observe lim^o ^^(1, c) = 2. (This is clearly correct 
since Rayleigh's equation becomes ip" = and the assumed boundary conditions imply ip = 1 + y, 
which when evaluated at y = 1 gives 2.) Defining F(c) := W(c)/[(U2 — c)(c — U\)\ gives F(c, 0) ~ 
2e~ 2z0 /R 2 , which is consistent with the A; = result of [5]. 

The principle of permanence can be demonstrated explicitly by using the solution to Rayleigh's 
equation written as ip(y,c) = Y^=o VVi (Ui c ) > where 

VWiG/,c) = /jy-y') (&+ ^l^ n {3/,c)dy'. (63) 

Using Cauchy's inequality, it is not difficult to prove that this series converges uniformly for all 
y ^ y c . If we choose ipQ so that the boundary conditions of (2) are satisfied, evaluate (63) at y = 1, 
and slip the c — > oo limit through the integral sign, then we obtain <J>^?(y). 

As an aside, note that each ^ n is analytic in c. Thus because of uniform convergence, ip(l,c) 
is analytic for c 7^ £7(1). This is true even for profiles U(y) that are not analytically continuable 
into the complex plane; it is only necessary for the integral in (63) to exist to get analyticity in 
c. This puts a relatively mild restriction on U . For example, if c is not in the flow domain, then 
yU"(y) G Li[-1, 1] is sufficient. 



The exterior legs: 3^4 and 8 — > 1. 

We show below that if c is not in the flow domain, which is the case on the legs 3 — ► 4 and 8 — > 1 
neither $<(y, c) nor <J>>(y, c) can vanish. Evaluating at y = 1 we see the same is true for W. 
We also include a proof that W is a monotonic functions of c on 3 — > 4 and 8 — > 1. Thus these 
pieces of the contour map into curves that cannot cross into the right hand portion of the W-plane. 

To prove the above statements we use a formula introduced by Green in the first half of the 
nineteenth century (see e.g. [22] or [10]). For Rayleigh's equation, Green's formula is 

which is valid for any solution ip. This formula can be derived from Rayleigh's equation by multi- 
plying by ip, manipulating, and integrating. It is important to remember that c is assumed to be 
real and outside the domain of integration. 



dy' > , (64) 
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Upon taking yo = — 1, -ip = <!><, and applying the boundary condition $<(— l,c) = 0, (64) 
implies 



>o, 



(65) 



U(y)-c 

for all y <G ( — 1, 1] and c [Lq, l^]- This inequality means that neither factor can vanish for y in the 
interior of the flow domain. We know that $><(— 1, c) = and that $><(— 1, c) = 1, and therefore in 
a neighborhood of y = —1, by continuity of the solution, 3><(y, c) > 0. Thus $<(y, c) > for all 
y <G (—1, 1], and inequality (65) implies that the [ ]-factor must also be positive. (Note that neither 
factor can be singular by the existence theorem applied to Rayleigh's equation.) We mention, for 
later use, that a similar argument shows $>(y,c) > 0. 

Evaluating <!>< at y = 1 yields W(c) = — $<(1, c) < for all c <G (— oo, Lq), which is our desired 
result for the leg 3^4. Similarly, W(c) < for all c G (U 2 ,oo), our desired result for the leg 
8-» 1. 

We now further demonstrate that W is monotonic, although this is not strictly needed for the 
proof. To this end we differentiate Rayleigh's equation with respect to c, yielding 



dip" 
~dc~ 



- k 2 + 



U"(y) \diP 



U"(y) 



■1>, 



U(y) - c) dc [U(y) - c] 2 
which by the method of variation of parameters, is seen to have the following solution: 
d^(y,c) [v V(y',c) U"(y') 



dc 



ipi(y',c)ip 2 (y,c) -^ 2 {y\c)^i{y,c) dy' . 



-1 W 12 [U(y')-cf 
Here ip\ and ip\ are any two independent solutions of Rayleigh's equation and 

W 12 (c) := tp! (y, c) ip' 2 (y, c) - ip 2 {y,c)ip[ {y,c) . 



(66) 



(67) 



(68) 



Observe that (67) satisfies dip(— l,c)/dc = 0, which is consistent with ip(— l,c) = 0, and thus we 
may assume ip(y,c) = <5<(y, c). Letting ip\{y,c) = <3?<(y,c) and ip 2 (y,c) = <£>(y, c) gives 



W 12 (c) := [^(y.cj^^c) - $>(y, cX(y, c)J = -$>(-l,c) = *<(l,c) 
and the expression (67) implies 



^<(l,c) = 
dc 

Integrating (70) by parts gives 

3*<(l,c) = | : 



C/"(y' 



^<(y / ,c)$ > (y / ,c)dy / . 



(69) 



(70) 



-2C/ ,: 



dc 



[U{y')-cf 



$<(y')$>(y') 



+ 



u> 



[^(yOM^O + MiO^O]} dy' , (7i) 



[t/(y')-c] 2 

which upon insertion of (69) into its second term can be manipulated into 

d*<(l,c) 



dc 



K( C )+/x(c)$<(l, C ), 



(72) 
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where 



and 



M (c) := £ 



K(c) 



[£%') - c 
[W)-c] 4 



■dj/' = 



<f<(y')- 



U2-U1 

(f/i-c)([/ 2 -c) 

W) *<(*/') 



(73) 



(74) 



Clearly for c ^ (U\,Uz), /x(c) > 0. The [ ]-f actor of (74) is precisely the [ ]-f actor of (65), which 
we showed is positive. Since U'(y) > and <3>>(y) > 0, for y G [—1, 1), we see that K(c) > 0. Thus 
we have established monotonicity: <?$<(!, c, k) / do or cW(c, k)/dc < for c ^ (£/i, L^)- 



The boundary regions: 4^5 and 7^8. 

Now we consider the pieces of the contour that skirt the boundaries of the flow domain. These are 
the pieces 4^5 and 7 — ► 8 of Fig. 2. We use the expressions i/j g and of (16) and (17) to construct 
<!>< (y, c) according to the procedure described in Section 4. Upon enforcing the boundary conditions 
as described in that section, it is not difficult to show that the Wronskian has the following form: 

W(c) = -$<(1, c) = Xi(c)(C/! - c) + X 2(c)(U 2 - c) 

+ (U 2 - c){U! - c)[ X3 (c) In(l7i - c) + X4 (c) ln(C/ 2 - c)] , (75) 

where the functions Xi for i = 1 — 4 are analytic functions of c. The piece 4^5 can be pulled down 
into a little semicircle on which c = U\ + 5e l6 ', where S > and 9 € (— 7r, 0), which is still consistent 
with c being in the upper half plane. Evidently, lim^o W(c) = X2(Ui)(U2 — U\). Since this piece 
of the contour must connect to the piece that terminates at 4 as 3 — > 4, it follows that the piece of 
the contour 4^5 maps into a point on the negative real axis of the VF-plane. Similarly, on 7 — > 8, 
where c = U2 + Se td , 5 > 0, and 6* £ (— ir, 0), we obtain lim^o^( c ) = Xi(^ r 2)(^ r i — c^), and the 
piece of the contour 7 — > 8 maps into a point on the negative real axis of the JU-plane. 

So in conclusion, since we have proven that for c £ C the real part of W is negative, the only 
part of the contour that can give rise to a change in the argument of W is the piece from 5^7. 
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Figure Captions: 



Figure 1: Nyquist plots for the single inflection point profile, U(y) = tanh/3y. (a) Four plots for 
(3 = 2 and k = 0, 1, and 2, and the critical value for the onset of instability, k = k c ~ 1.832 
(dashed curve), (b) Four plots for k = and (3 = 1, 1.5, and 2, and the critical value for the 
onset of instability, (3 = (3 C ~ 1.2 (dashed curve). 

Figure 2: The contours C, C , and -ff in the c-plane. The closed contour C runs the entire circuit 
from 1— >2— »8— >1, with the portion from 3 — > 1 lifted infinitesimally above the real 
axis. The contour C is C with the piece along the flow domain, 5^7, removed. The contour 
H (dashed) is Howard's semicircle, within which the unstable eigenvalues must lie. 

Figure 3: Stability boundaries for the three equilibrium profiles of Sec. 7. Figures (a) and (b) 
depict the nucleation of instability through k = 0, while (c) depicts the nucleation of instability 
through k ^ 0. In (b), the stability boundary of the second and third unstable modes is also 
shown. 

Figure 4: Nyquist plots for the multiple inflection point profile, U(y) = y + (sin7y)/27, with 
k = 1. (a) The case 7 = 5 has one encircling of the origin and thus one unstable eigenvalue, 
(b) In the case with 7 = 7.5 two more loops have appeared and the curve almost encircles the 
origin two more times. That is, it is close to the threshold of the instability of two further 
modes, (c) The case 7 = 10 has three encirclings of the origin and three unstable eigenvalues, 
(d) Depiction of the monotonic dependence of the crossing values to the left of the asymptote 
(1,0) as a function of k. 

Figure 5: Vorticity and vorticity gradient profiles corresponding to the velocity profile U(y) = 
y + 5y 3 + /tanh4(y — 1/2) with equilibrium parameter values / = 1.3, / = 1.65, and / = 2. 
At / = f c « 1.59 two inflection points emerge. 

Figure 6: (a) Nyquist plots for the profile U(y) = y + 5y 3 + /tanh4(y — 1/2) for k = 0, 0.75, and 
1.5, and / = 1.62. (b) Magnification that indicates the onset of instability through k / 0. 

Figure 7: Plots of e r (cj) against k at / = 1.62, where c/, / = 1,2, are the two inflection points 
leading to crossings of the Nyquist plot in Fig. 6 to the left of the asymptote (1,0). 
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